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The basic ingredients of osmotic pressure are a solvent fluid with a soluble molecular species which 
is restricted to a chamber by a boundary which is permeable to the solvent fluid but impermeable 
to the solute molecules. For macroscopic systems at equilibrium, the osmotic pressure is given by 
the classical van't Hoff Law, which states that the pressure is proportional to the product of the 
temperature and the difference of the solute concentrations inside and outside the chamber. For 
microscopic systems the diameter of the chamber may be comparable to the length-scale associated 
with the solute- wall interactions or solute molecular interactions. In each of these cases, the assump- 
tions underlying the classical van't Hoff Law may no longer hold. In this paper we develop a general 
theoretical framework which captures corrections to the classical theory for the osmotic pressure 
under more general relationships between the size of the chamber and the interaction length scales. 
We also show that notions of osmotic pressure based on the hydrostatic pressure of the fluid and the 
mechanical pressure on the bounding walls of the chamber must be distinguished for microscopic 
systems. To demonstrate how the theoretical framework can be applied, numerical results are pre- 
sented for the osmotic pressure associated with a polymer of A'' monomers conflned in a spherical 
chamber as the bond strength is varied. 

PACS numbers: 05.20.-y,05.70.-a,05.10.-a 



I. INTRODUCTION 

Osmotic effects are thought to play an important role 
in many physical systems associated with biology and 
technological applications. Examples in technological ap- 
plications include microscopic devices designed to pump 
fluids P, [2, [1| , actuate forces through swelling 0, Q , or 
deliver drug doses [1, 0j 9 ■ Some macroscopic osmotic 
mechanisms in biology include the exchange of blood con- 
stituents in capillaries with surrounding tissues ^9] and 
the processing of fluids in tissues of layered epithelial cells 
in the intestines and kidneys [l^, [ll| . 

At a more microscopic level, individual cells contain 
a high concentration of charged molecules, in which os- 
motic effects must be controlled actively to avoid exces- 
sive swelling and bursting of cellular structures [l^ . 
In fact some of the mechanisms by which neural cells 
transmit electrical signals through the production and 
propagation of action potentials may have as their evolu- 
tionary origins the pumping mechanisms developed by 
cells to use counter-ion fluxes to compensate for the 
harmful effects of osmotic pressure [l^, [ill- Other ex- 
amples include the study of the pressures involved in 
DNA confinement in virus capsids [l5[ , packaging of pro- 
teins in small cellular vesicles in cell organelles [lal, and 
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even mechanisms of gel swelling in propulsion in micro- 
organisms such as myxobacteria [3, 13, [l3| ■ 

With new experimental techniques, such as optical 
trapping and molecular tagging, it is now feasible to ob- 
serve and measure forces and displacements in systems 
on the length-scale of hundreds to tens of nanometers [H, 
ITqI . [21I I , and the underlying physical processes can begin 
to be explored quantitatively at very small length-scales 
in biological and synthetic systems [isl. [isl. [Tol. [20. [2^ [23t . 
The modeling and analysis of these types of systems mo- 
tivates using a theoretical framework from statistical me- 
chanics which is general enough to apply to these systems 
without the usual type of hard-wall and other scale sep- 
aration assumptions appropriate in classical thermody- 
namic systems. 

In this paper we shall discuss a general microscopic the- 
ory for osmotic pressure at equilibrium and draw some 
contrasts with the classical theory of van't Hoff [2^ . 
First, in Section HIl we introduce two statistical quanti- 
ties to characterize pressure associated with osmotic phe- 
nomena: one related to the pressure built up in the fluid 
and another to the pressure felt by the confining wall of 
the chamber. We find in Section [1111 that while both lead 
macroscopically to equivalent notions of osmotic pres- 
sure, they differ in general for microscopic systems with 
interaction length scales comparable to the chamber size. 
Each form of osmotic pressure has theoretical connections 
to macroscopic thermodynamics; the osmotic wall pres- 
sure maintains its statistical mechanical relationship to 
changes in free energy with respect to volume changes, 
while the osmotic fluid pressure is directly related to the 
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notion of osmotic pressure in macroscale nonequilibrium 
thermodynamics. A specific example illustrating the dis- 
tinctions between the two notions of osmotic pressure 
and their deviations from the classical van't Hoff law is 
presented in Section llll Bl and Sect ion [TVl in which nonin- 
teracting solutes have a solute-wall interaction potential 
given by a power-law of the distance of the solute to the 
wall. 

We then proceed in Section|V]with an example with so- 
lute particles bound by "string" forces in a hard-walled 
potential to illustrate how the interaction length scale 
between solute particles affects the osmotic pressure. Af- 
ter preparing in Section IVII with an alternative expres- 
sion for the osmotic pressure in terms of corrections from 
the van't Hoff law, we illustrate in Section IVIII the the- 
oretical formalism through numerical Monte Carlo cal- 
culations for the osmotic pressure for a polymer con- 
fined within a spherical chamber with various binding 
strengths (and therefore bond lengths). The theoreti- 
cal framework presented in this paper is expected to be 
applicable to model osmotic effects in many microscopic 
systems at equilibrium, and provides a step toward de- 
veloping a non-equilibrium (or near-equilibrium) theory 
for microscopic osmotic phenomena. 



II. MICROSCOPIC OSMOTIC PRESSURE 

We begin by posing two possible means for describing 
the pressures associated with osmotic phenomena at a 
microscopic level in a chamber with a confining boundary 
(wall) permeable to the solvent but not the solute. We 
are concerned here with developing a microscopic anal- 
ysis, rather than a macroscale thermodynamic descrip- 
tion. In particular, we are interested in a theory based on 
quantities that could be measured in a microscale exper- 
iment or computed from a microscale simulation, using 
a numerical method for exarnple such as the stochastic 
immersed boundary method [25| or Stokesian dynamics 

The first quantity we shall consider is the mechanical 
pressure exerted by the solute on the wall. The second 
quantity is the hydrostatic pressure built up in the sol- 
vent in the interior of the chamber. While these two 
notions lead to the same pressure quantities for macro- 
scopic systems, referred to as the osmotic pressure, we 
shall show that different results are obtained for the pres- 
sures when considering microscopic systems. 

The osmotic pressure definition based on the wall pres- 
sure maintains some classical thermodynamic relation- 
ships involving pressure and other statistical mechanical 
quantities in nonideal systems for which the length scale 
of the wall and solute interactions are nonzero and fi- 
nite. On the other hand, the osmotic pressure definition 
based on the manifested fluid pressure can be shown to be 
equivalent, in a certain local sense, to the osmotic pres- 
sure definition used in the theory of macroscopic nonequi- 
librium thermodynamics [27l [2g | to describe the driving 



force behind solute fluxes across a semipermeable wall or 
membrane. 

Both notions of osmotic pressure therefore have some 
underlying theoretical justification in both microscopic 
terms and in relation to macroscopic thermodynamics, 
though we will show through theory and example that 
they are not equivalent. Both types of pressure moreover 
would appear to have practical relevance for modeling 
and analysis depending on the application. 

When investigating swelling phenomena and the elastic 
forces of a confining membrane, for example, a notion of 
pressure involving the average force exerted per unit area 
on the confining wall may be of particular interest. In 
contrast, when investigating the role of fluid dynamics 
in transport and dilution of the concentration of solute, 
for example, the hydrodynamic pressures which drive the 
flows may be of more relevance. 

In what follows, we generally restrict attention to the 
case of N identical and possibly interacting solute par- 
ticles confined in the chamber. Furthermore, we shall 
assume throughout that the potential energy of the par- 
ticles is of the form: 

JV 

*(X1,--- ,Xjv) = F(X1,--- ,XAr)+^$(Xfc). (1) 

fc=l 

In the notation, V{x.i, • • • , xjv) models the energy of the 
solute-solute interactions and ^(x) models the interac- 
tion of a solute particle with the chamber wall. The Xfc 
denote the spatial coordinates for the fcth solute parti- 
cle. Much of the theory can be readily extended to more 
general potentials. 



A. Osmotic Wall Pressure 

The first notion of "osmotic pressure" we consider will 
be defined in terms of the average forces that solute par- 
ticles exert on the walls of a confining chamber. We shall 
assume that the solute- wall interaction potential $ arises 
from a uniform areal distribution of particles on the wall 
which have an isotropic interaction with the solute, given 
by a potential 4>. More specifically, we shall assume that 
for a confined solute particle at interior location x the in- 
teraction force exerted on the differential area dy of the 
bounding wall can be expressed as: 

G(x-y)rfy (2) 

where G(r) — — Vr(/)(|r|) and (f> is the isotropic interac- 
tion potential of the wall particles with the solute parti- 
cles. 

The force acting on a solute particle in the chamber 
interior which arises from the wall interactions is then 
given by: 

Fwi(x) = / -G(x-y)dy (3) 
Jdn 
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where J7 defines the container volume and the con- 
tainer boundary, which we often simply call "a wall." We 
discuss in Appendix |A] how the form of the solute-wall 
interaction (G or (/)) can be inferred from an observed or 
known structure for the wall force F^i or potential $. 

To keep focus on the main issues of interest, we will 
treat the container shape as prescribed and fixed. The 
results would apply to more general systems with a flex- 
ible boundary (such as a membrane-bound vesicle), pro- 
vided the thermal fluctuations of the bounding surface 
in equilibrium could be neglected due to their amplitude 
or time scale. In this case, we would not be predicting 
self-consistently the shape of the container, but simply 
taking its (time-averaged) shape in thermal equilibrium 
as an input to deflne fl. 

Forces of the form Eq. ([3]) can arise in a system in a 
variety of ways. For example, a solute particle (possibly 
a macro-ion) in the interior of a chamber could inter- 
act with like-charged particles composing the wall sur- 
face [2^. In this case, G would represent the screened 
Coloumbic force per unit area exerted on the surface as a 
consequence of the charge density of the wall and solute 
while F„i would be the screened Coloumbic force acting 
on the solute particle when it occupies location x. In 
principle, both of these forces could be derived from a 
common electrostatic potential for the system 'si']. An- 
other example would be for G to represent an effective 
force of interaction between a particle and the wall aris- 
ing from time-averaged steric interaction forces such as 
those between a particle and a polymer brush coating 
the wall [131 ■ One could even consider within this frame- 
work the degrees of freedom of a single monomer of a 
large, compactly folded polymer as constituting an effec- 
tive particle confined by the time-averaged effects of the 
rest of the monomers in the polymer. 

To define a pressure for these types of systems we can 
consider the average forces which act in the normal di- 
rection across the wall surface: 



Pwl = 



/ 1 



\\dn\ 



I f]G(Xbl(t)-y).nyrfy 

Jdn .J / 



(4) 



where \dVl\ denotes the surface area of bounding wall 957, 
fly denotes the unit outward normal to the surface at y 
and (•) denotes a time average over the position of the 
N solute particles {X[^l(t)}|^^, which we shall assume is 
equal to the statistical ensemble average. In the notation, 
Xt-'l [t) denotes the position of solute particle j at time 
t. We shall refer to this definition of pressure as the 
"osmotic wall pressure." 

We can alternatively express this osmotic wall pres- 
sure in terms of the density of the wall-solute interaction 



potential as: 

/ 1 



/ 1 



J Oil Jn 



\\dn\ 

•(5(x-X[^l(t))rfxdyJ 

where we have introduced the Dirac delta function 5(x) 
in order to express the statistical average in terms of the 
solute concentration |32l |: 



/ 

c(x) :=(^5(x-XW(t)] 



(6) 



This concentration is proportional to the probability den- 
sity for any of the particle positions X'^l (t) to be at lo- 
cation X when the system is in statistical equilibrium. 



B. Osmotic Fluid Pressure 

Another approach to studying osmotic pressure is to 
consider the solvent fluid and, in particular, how the 
forces of interaction between the solute and wall are man- 
ifested in the hydrostatic pressure. For a Stokesian fluid 
we have that the local fluid velocity u satisfies: 



du 



V-u = 



= -/iAu - Vp f ft 



th 



(7) 
(8) 



where p is the fluid density, jjL is the dynamic viscosity, 
p is the fluid pressure, f is the force density exerted by 
the solute particles, ft-h is the force density associated to 
thermal fluctuations [2^, [H, [sH] • 

Working with the immersed boundary method approx- 
imation 36] for the interaction of a fluid with N identical 
solute particles, we express: 



N 



f(x,i) = E 



-Vj^'(AJ')(5(x-XW(i) 



(9) 



In the notation, 'X.^\t) denotes the position of solute 
particle j at time t and V^- denotes a gradient with 
respect to the spatial coordinate of solute particle j. For 
compactness in the notation we also define convenient 
expressions involving composite vectors of the particle 
positions as follows: A^ = AJ(x, {XW(<)}f^i_fe_^j-) = 



[xW,...,X[-'-il,x,X[^'+il 
A,(x,{x,}f^,_,^^.) 



,X[^1] 



and Ao 



I Xj — 1 , X, Xj_|-i , 



The (5(x) describes the manner in which the force on 
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the solute particles is distributed to the fluid; for small 
idealized point particles this would just be the usual 
Dirac ^-function. 

More general particle-fluid interaction functions, how- 
ever, could be used both for modeling and numerical rea- 
sons, as in the immersed boundary method [s^. Our end 
results will generally be posed in a manner which can be 
appropriate for 6 functions with zero or finite width, but 
for conciseness we will treat the 6 functions in deriva- 
tions as classical Dirac delta functions. Moreover, one 
could more rigorously describe the fluid-particle interac- 
tions in terms of rigid or flexible moving boundaries, but 
we choose to use the immersed boundary approximation 
with Dirac S functions in this initial work since it will con- 
vey our central ideas and distinctions with a minimum of 
technical distraction. 

Since no external driving force is applied to the cham- 
ber we have at statistical equilibrium that the average 
fluid velocity must be (u(x)) — 0, [sj, which, along with 
the fact that the thermal force density fth has zero mean, 
implies 



= - (Vp) + (f> 



(10) 



Expressing the force density f explicitly in terms of the 
potentials of the forces acting on the solute particles, we 
can define an average pressure gradient: 

Vp(x) (Vp(x)) - (f(x)) 

/ N 



- (5]-Vx<i>(x)<5(x-X[-'"l(<) 



/ ^ 

+ /^-VjT/(AJ>(x-Xb'l(t)^ 

= -V,,$(x)c(x) 
/ ^ 

+ /^-VjF(AJ)(5(x-X[-'l(t)^ 

„ N \ 

n (Hx.-xi'i«)).x.)\ 

= -Vx«'(x)c(x) 

W „ N 

+ E / -Vx,y(A,)p^(A,) n 

(11) 

In the notation, c(x) is the local concentration density of 
the solute as defined in Eq. ^ and 

p^(xi, X2, . . . , xat) = (l[ 5(x, - Xbl (t))^ (12) 

is the N particle correlation function in thermal equilib- 
rium. All terms in the last sum are equal if the particles 
are truly identical and exchangeable, but we leave the ex- 
pression as an explicit sum to incorporate systems with 



bonded interactions. For example, in Section FVllI we con- 
sider a confined polymer in which the monomers on the 
ends have a different bonding structure than monomers 
in the interior of the polymer chain (ends have only one 
bond, interior monomers have two). 

By integrating the pressure gradient along a path start- 
ing from a reference point Xyi outside the domain 
(where we set the pressure to the reference value of zero), 
we can thereby define what we will call the "osmotic fluid 
pressure" Pfl(x) inside the domain: 



Pft(x) := p(x) = / Vp(x) • dx 



(13) 



To compare this notion of pressure with the osmotic wall 
pressure we shall mostly present the osmotic fluid pres- 
sure at the center of the domain (which we shall take to 
be x = 0), but the definition of the osmotic fluid pressure 
field extends throughout space. 



III. STATISTICAL MECHANICAL FORMULAS 
FOR OSMOTICALLY RELATED PRESSURES 

As a starting point for our ensuing analy- 
sis, we show how both osmotically related pres- 
sures Eq. (HI) and Eq. defined in Section |TT] can be 
represented in terms of the partition function for the 
solute particles in a soft-walled potential. In Subsec- 
tion [1112 we will verify that for noninteracting particles 
in the hard-wall limit, we recover the standard van't 
Hoff's law for osmotic pressure both on the wall and in 
the fiuid. 

We consider first the pressure Eq. ^ induced by the 
solute against the chamber wall. In what follows, we will 
relate the osmotic pressure to changes in the Hclmholtz 
free energy (or equivalently the partition function in the 
canonical ensemble) under the deformation of the con- 
tainer volume, and the formulas will only take familiar 
statistical mechanical form if we scale the wall-particle 
interaction force (and energy) density inversely with the 
local area of the surface. This corresponds to the physical 
situation in which the number of particles making up the 
wall, rather than their areal density, remains fixed under 
deformation. Therefore, we define a one-parameter fam- 
ily of deformations (diffeomorphisms) rj^l-) which deform 
the chamber boundary dQ = dCl{sQ) into a new bound- 
ary dfl{s) by extending or contracting the surface at a 
constant rate along the local outward normal vector liy 
at each y G dn{s). The differential equations describing 
this family of deformations are: 



dvsiy) 



ds 



(14) 
(15) 



These mappings are generally smooth only over a finite 
interval of parameters s containing sq. Then, taking 
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our (f) and G to refer to the potentials and forces as- 
sociated with the reference chamber surface 90(so), we 
define the potential for solute-wall interaction forces in 
the deformed chamber boundaries through: 

$(x,s) = / 0(|x-y|)det(V777'(y))dy, (16) 
Jdn(s) 

where the Jacobian factor det(Vj77^(y)) is inversely pro- 
portional to the local expansion of area under the de- 
formation. The potential so defined is naturally related 
to the force defined above in Eq. ^ through F(x) = 
— Vx$(x, So)- The total normal force acting outward on 
the surface of the wall exerted by a solute molecule at 
location x is then given by: 



/i(x) = / -[ny.Vy0(|x-y|)]det(Vr7:i(y))dy 

Jdn(s) 

[ ny • Vx0(|x-y|)det(Vj7;i(y))dy 

( n^^(y.) • Vx0(|x-77,(y')|) dy' 
dVsiy') 



ds 



Vx<^(|x-r,,(y')|)dy' 



I- / 0(|x-r,,(y')|)dy' 

Jan{so) 



c»$(x,s) 
ds 



(b{\ic-y\)det{Vr]-\y))dy 



(17) 



The solute concentration c(x) is obtained as a one- 
particle contraction of the n-particle correlation function: 

c(x)=^ / pn{A,) Yl dxfc, (18) 

which, in thermal equilibrium, can be expressed in terms 
of the Boltzmann distribution: 



pjv(xi,X2, • • • ,XAr) 

= Z(s)"^exp 



1p{-Kl,X2, ■ ■ ■ ,Xjv) 



(19) 



where T is the temperature, is Boltzmann's constant, 
and the partition function for the solute particles con- 
fined by a surface d^{s) is given by: 



Z{s) — I exp 



exp 



V^(xi,X2, ■ ■ ■ ,Xjv) 
y(xi,X2, ■ ■ ■ ,XAr) 

knT 



N 



fe=i 



exp I , ^ J ^ ' j Yi ^^fe- 

/ k=l 



knT 



(20) 



Substituting Eq. HH]) and Eq. ^ into Eq. ©, we obtain 
for the osmotic wall pressure: 



p 



^|^Mx)c(x)dx 

1 /■ 9$(x,s) 
W\Jn ' 



N 



•P7V(X1,X2, • • • ,Xjv) Y[ ^^fc 

fcsT d 



\dfl\ ds 



fc=i 

ln(^(s))L=.„. 



(21) 



and in the last expression, we have indicated that the 
pressure is to be calculated for the initially specified 
chamber wall s = Sq. Our result Eq. (pij) agrees with the 
thermodynamic definition of pressure as the (functional) 
derivative of the Helmholtz free energy —ksThiZ with 
respect to volume V along our one-parameter family of 
deformations, for which dV = |9ri|ds. We remark that 
since the bounding walls are assumed to be impermeable 
to the solute, the confinement energy must diverge at 
the boundary: $(x, s) ^ oo as x ^ 9ri(s), which has 
the effect of restricting integration to only the interior 
of the chamber. The reason that the osmotic wall pres- 
sure is only equal to the volume derivative of the free 
energy along a certain one-parameter family of shape de- 
formations appears to be related to the fact that nonlo- 
cal interactions (such as by a soft wall potential) make 
the local pressure no longer constant within a container, 
and therefore along the surface, so that different kinds of 
shape deformations will have different free energy costs 
per volume of deformation. 

Next we develop a statistical mechanical representa- 
tion for the osmotic fluid pressure Eq. p^ . our sec- 
ond notion of osmotic pressure. We shall again assume 
that the system at statistical equilibrium has Boltzmann 
statistics Eq. (fTOl) , from which it follows that the expres- 
sion for Vp(x) in Eq. (jlip can be simplified as follows: 



N 



Vp(x) 



ksT 
Z 

ksTVcix). 



exp 



^(A.) 
knT 



N 



k=i,k^j 



From Eq. ^ and Eq. 



we obtain 



Pfl(x) = A:bTc(x) 



(22) 



(23) 



This shows that our form of the osmotic fluid pres- 
sure connects with macroscopic notions in nonequilib- 
rium thermodynamics [27j . 
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A. Hard- Wall Limit 

We now show how these formulas relate to the van't 
HofF law for osmotic pressure both on the confining wall 
and in the fluid in the case that the particles are nonin- 
teracting and the confining potential has a classical hard- 
wall potential of the form: 

0, X e n{s) 

OO, X ^ f2(s) 



$(x, s) = 



(24) 



Substitution of this "hard-walled" confining potential 
into Eq. ((20|) . we have that Z{s) is just the iVth power 
of volume V{s) of the chamber f2(s). By Eq. (PT|) . the 
osmotic wall pressure is then: 



(25) 



where Cg = N /V{s{)) is the concentration, recovering the 
well-known van't Hoff's Law equation 24]. Similarly, 
from Eq. (|23p . we have: 



ft(x) 



NkBT{V(.so)) 



NkeT 



and again find that the classical van't Hoff's Law is re- 
covered in the fluid pressure. 

We remark that the key to obtaining the classical van't 
Hoff's Law in both cases was to consider the hard-walled 
limit of the confining potentials which arises naturally 
when the length scale of particlc-wall interactions are 
very small relative to the diameter of the confining cham- 
ber. In this regime a theory similar to our approach was 
developed in [J] , in which a careful analysis is made of the 
balance of mechanical forces arising from the solute inter- 
actions with the walls, again under the assumption that 
the chamber diameter is much larger than the length- 
scale of the particle interaction forces. While this as- 
sumption typically holds for macroscopic systems, when 
the chamber size becomes sufficiently small, the results 
of this limit are no longer strictly valid and corrections 
are required to the classical theory reflecting some of the 
microscopic features of the system. 



B. Steric Interactions with the Chamber Walls of 
Non-negligible Length-Scale 

For microscopic systems, the length-scale on which 
the solute particles interact with the wall may be non- 
negligible relative to to the diameter of the chamber. For 
example, if a rigid spherical particle of radius £ is con- 
fined in a spherical chamber having radius R, the steric 
interactions would confine the particle center to have only 
configurations with |x| < R — I. We shall now illustrate 
how this effects the osmotic pressure when £ is compara- 
ble in size to R. 

More precisely, in the case of a spherical chamber of 
radius R the confinement potential is of the general form: 



<i>(x,i?) = / 
■Id 



|y-x|)||dy. 



(27) 



where (j){p) denotes the solute-wall potential for a given 
reference sphere of radius Rq and det(V^7^(y)) = 
Now for R^ ^ R- £ < Rwe have: 



0, 

oo, 



p>£ 



(28) 



where p = |y — x|. By radial symmetry this gives for 
r — |x| the solute confinement potential: 



$(r, R) 



0, r<R~£ 
oo, r > R- £ 



(29) 



The partition function and its derivative are obtained 
explicitly as: 



Z{R) = -{R-ef 



dZ{R) 
dR 



= A7r{R-£)^ 



(30) 
(31) 



fcsTco, (26) From Eq. (fT5|) . the osmotic fiuid pressure is given by: 



ft(0) = 



f (i? 



£Y 



(32) 



From Eq. (j2ip . the osmotic- wall pressure is given by: 



P. 



wl 



Z 



knT 



4f (i? 



£f 



(33) 



where £ < R. 

The difference in these two notions of pressure in this 
case can be interpreted geometrically. In particular, the 
confinement forces restrict the solute particles within a 
spherical region of radius i?^ = R ~ £. This has the 
effect of generating within the confinement region the 
same local confinement forces as a hard-walled poten- 
tial with a wall occupying the spherical shell of radius 
Ri, — R — £ and yields the same osmotic-fluid pressure 
as in the hard-walled case. However, the particle-wall in- 
teractions, which generate the confinement forces, occur 
from a wall occupying the spherical shell of radius i?, as 
opposed to a spherical shell of radius i?* as would occur 
in the strictly "hard- wall" case. From the principle of 
equal and opposite forces, the particle-wall interactions 
now exert forces in the normal direction over a wall with 
greater surface area than in the short-range hard- walled 
case, thus reducing the osmotic wall pressure by the ge- 
ometric factor (1 — £/ R)^ = (i?*/i?)^, which is the ratio 
between the surface areas. 

While this example is rather special, it illustrates 
clearly one mechanism by which differences with the clas- 
sical theory of osmotic pressure can arise. This also il- 
lustrates how notions of osmotic pressure in terms of the 
wall pressure and fluid pressure can differ markedly when 
the length scales of the particle- wall interactions become 
comparable to the diameter of the chamber. We next 
demonstrate how osmotic pressures behave for smooth 
potentials having long-range particle-wall interactions. 
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IV. OSMOTIC PRESSURE FOR 
NON-INTERACTING SOLUTE PARTICLES 
CONFINED BY A SOFT- WALLED POTENTIAL 
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We now examine the behavior of the osmotic pressure 
for a system of non-interacting solute particles which in- 
teract with the chamber wall through a smooth long- 
range potential. For simplicity, we will consider the os- 
motic pressure exerted by a single solute particle; the 
case of N non-interacting solute particles of course sim- 
ply multiplies the pressures by a factor of N. We shall 
consider here the class of repulsive potentials for a spher- 
ical chamber of radius Rq of the form: 



0(p) = Cp-", 



(34) 



which, under volume dilations to new radii R, induce 
effective wall potentials 



$(r, R) 



We consider only a > 2 to ensure that lim^^ij $(r, R) 



(35) 



The partition function and its derivative can be ex- 
pressed as: 



Z{R) = inR^ / exp - 



-(1 



(a - 2)Ar 



r^dr 



(36) 



and 
dZ_ 
dR 



3Z{R) 

R 
47ri?2 

A 



a 



a - 2 
exp 



1 /(I -r)2-" - (1 +r)2 



(1 - r)2-" - (1 + r)' 
(a-2)Ar 



r^dr 



(37) 



where A = ksTR" /2ttCRq can be regarded as a reduced 
temperature. These expressions were obtained by us- 
ing Eq. (flB]) and Eq. ([20|) with R — s and making the 
change of variable r = Rr' . This gives the osmotic fluid 
pressure: 



^fl(O) 



Qi(A,«) 



and the osmotic wall pressure: 



Pwl — 



Q2(A,a) 



(38) 



(39) 



where 



)i(A,a) = exp ( -- 



3 / exp 
Jo 



(l-r)2 



(l + r)2 



{a - 2)Xr 



r^dr 



(40) 



10 



■a 10" 



10" 



O 



10" 



10 



10 r 



a = 3 



a = 5 



a = 40 



10 



10 ■ 



10" 



^/x^' 



FIG. 1: Qi Correction Factor for Soft- Walled Potential. For 
A oo we have and Qi 1. This corresponds to 



_ the plotted horizontal lines. For A ^ we have 1/A 



1/2 



and Qi/{{a - r 
the plotted diagonal line on the right. 



^\^/'^^ . This corresponds to 



and 



2(A,a) = 1 + 



a-2 



(41) 



1 / (l_^)2-._(l + ^)2-. \ 

Jo *^^P (a-2)\r ) 



•■^dr 



For this system with the power-law potential, we can 
express Qi compactly as 



(42) 



and as 



Q2(A,a) = 1 



1 



a 
3 



3kBT 



(l-r) 



2-a 



(1 



\2-a 



(a - 2)Ar 



(43) 



The last expression shows that Q2 is one plus a term 
proportional to the ensemble average of the energy of 
the system over all configurations of the confined particle 
with respect to the Boltzmann distribution. 

We observe first of all that the soft-walled nature of 
the potential produces corrections to the van't Hoff law 
for the osmotic pressure of the chamber in both the os- 
motic wall pressure and fluid pressure. To compare these 
correction factors with each other, we first examine their 
asymptotic behaviors in the limits A — > and A — > cx). 
Note that the term A represents the ratio of the energy 
scale of the thermal fluctuations relative to that of the 
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10" 



MX" 



10 



10 



FIG. 2: Q2 Correction Factor for Soft- Walled Potential. For 
A cx) we have and Q2 1. This corresponds to 

the plotted horizontal lines. For A we have 1/A^^^ 00 
and Q-ijct ^ 2/3 A. This corresponds to the plotted diagonal 
line on the right. 



confining potential, and therefore these limits correspond 
to low temperature and high temperature limits, respec- 
tively. As A — > 00 we have Qi,Q2 — > 1, which recovers 
the hard-walled limit. On the other hand, for small A, 

3/2N 



9V37r 



\ia-l)aj 



and Q2 ~ 2q;/3A, show- 



ing a behavior similar to a system at low temperature 
in which the structure of the long-range particle-wall in- 
teractions plays a significant role. We note that the cor- 
rection factors for the osmotic fluid pressure and osmotic 
wall pressure diverge at different rates with respect to 
A as A ^ 0, indicating that the osmotic fluid and wall 
pressures will deviate significantly both from each other 
and from the van't Hoff law for soft-walled confining po- 
tentials with energy scales large compared to the temper- 
ature. On the other hand, plots of the behavior of the 
correction factors as a function of 1/A^/^ (see Figures [T] 
and [5]) show qualitatively similar behavior, in that the 
ratio of the osmotic pressures to the values given by the 
van't Hoff law increases monotonically from 1 as 1/A^/^ 
increases and diverges as A — > 0. 



V. OSMOTIC PRESSURE OF SOLUTE 
DIMERS CONNECTED BY "STRINGS" AND 
CONFINED BY A HARD- WALLED POTENTIAL 

We shall now discuss how the length scale of interaction 
between the solute particles can affect the osmotic pres- 
sures when this length scale is comparable to the cham- 
ber size. To illustrate through analytical formulas the 
basic mechanism by which this occurs, we shall consider 
the rather special case of two solute particles which are 



FIG. 3: Geometry of the Intersection of Two Spheres. The 
surface of any two spheres intersect in a circle carried in a 
common plane (dashed circle). We shall refer to this as the 
"plane of intersection" and further assume that the center of 
the smaller sphere lies within the volume of the larger sphere. 
For any such plane we may define three volumes by "slicing" 
the small and large sphere along this plane. We let ci denote 
the volume the sliced part of the larger sphere lying to the left 
of the "plane of intersection" . We let C2 denote the volume 
associated with the sliced part of the small sphere lying to 
the right of the "plane of intersection" . We let C3 denote the 
volume associated with the sliced part of the larger sphere ly- 
ing to the right of the "plane of intersection." To characterize 
the numerical values of the volume of these regions we define 
hi to be height of the spherical cap C3 and /i2 the height of 
the spherical cap C2. 



connected by a "string" and confined by a "hard-walled" 
spherical potential. More precisely, the potential energy 
of the system is given by 

*(xi, X2, i?) = $(X1, i?) + $(X2, i?) + l/(xi, X2) (44) 



where $ is the hard- wall potential given in Eq. . The 
potential V models the interaction between the two so- 
lute particles, located at Xi and X2, and is given by: 

T^/ \ f 0, if|x2— xi|<^ .... 
y(xi,X2) = <^ t-i! ■ ■ (45) 

^ ^ [00, otherwise ^ ' 

This potential can be given the physical interpretation 
of no coupling between the two particles until they at- 
tempt to separate past a distance £, at which point an 
infinitely strong restoring force constrains the particles to 
be a distance less than or equal £. This form of coupling 
is analogous to connecting the two particles by an inelas- 
tic "piece of string" of length £. The osmotic pressure for 
M noninteracting strings would of course just multiply 
the 1-string pressures we calculate by M. 
Using this form of the potentials we have: 



*(xi,X2,i?) = 



0, (xi,x2)er(i?,£) 
«), (xi,x2) ^r(i?,£) 



(46) 
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where 

r(i?,£) = {(xi,X2)||xi| < B, |X2| < R, |X2 -xil < £} 

(47) 

with r the admissible region in x R'^ in which both of 
the sohite particles are a distance no more than £ apart 
and contained within the spherical chamber of radius R. 
In this case, the partition function can be expressed in 
terms of a few quantities having a straightforward geo- 
metric interpretation and computed exactly. In particu- 
lar, the partition function 



Z{R) = / cxp 



^(xi,X2,fl) 



dxi dx2 (48) 



is the volume of the region r(i?, £) C R^. 

To compute this volume, we shall find it convenient to 
split the configuration space into two parts: T^^^ {R, £) = 
{|xi| < R- £} and T'^^^R,£) = {R - I < |xi| < i?}, 
where it is to be understood that r'^^^F'^^) C F. In the 
first region, for each xi the second particle X2 is free to 
assume values within the entire sphere of radius £ about 
Xi having volume ^t^. In the second region, for each xi 
the second particle X2 must lie within the intersection of 
the sphere of radius £ centered at xi and the sphere of ra- 
dius R centered at 0. The volume of the region so defined 
will be denoted by W*(|xi|). Using this decomposition, 
the partition function can be expressed as: 



Z{R) 



r(i) o 



r(2) 



Vl^*(|xi|)dxi 



47r 



An 



{R-ty 



(49) 



+ in Wi,{r)r^dr 
where the volume Wi, can be expressed as: 

W^r) = C{R,h,ir)) + —£^ ~ Cii,h2{r)) (50) 

with C denoting the volume of a spherical cap of height 
h for a sphere of radius r (see figure [S]) : 



The osmotic pressures are then given for ^ < i? by: 



^'fl(O) 



Pwl = 



2kBT 



. 16 \R) ^ 32 \R) 



. ^ 16 \R) ^ 32 \R) 



(54) 
(55) 



This gives the correction to the classical theory when the 
interactions of the solute particles becomes non-negligible 
relative to the diameter of the chamber. For £ > 2R, the 
particle interaction becomes trivial and both osmotic wall 
and fluid pressure assume the van't Hoff [l^l law value of 
2kRT/^TTR^. 



The above expressions verify that in the limit £ 



0, 



the osmotic wall pressure converges to the classical van't 
Hoff pressure ksT/ ^ttR^ of a single confined particle. 
This is to be expected as the two-particle string should 
coalesce into a single entity in the tight-binding limit. 
A similar transition between van't Hoff law behavior for 
the osmotic wall pressure exerted by dimers in the weak- 
binding and tight-binding limits was shown for spring- 
like dimers in numerical simulations with the stochastic 
immersed boundary method in (25j . 

The osmotic fluid pressure, by contrast, approaches 
the value 2kBT / ^ttR^ as ^ ^ 0, which corresponds to 
the van't Hoff law for two particles in the sphere, even 
though the two particles are tightly bound. This provides 
another example for how the pressure built up in the ffuid 
near the center need not be directly linked in magnitude 
to the pressure exerted by the solute particles on the 
chamber wall. 



VI. CORRECTIONS TO VAN'T HOFF LAW 
FOR INTERACTING PARTICLES IN 
SPHERICAL CHAMBERS WITH SOFT- WALLED 
POTENTIALS 



C{r,h) 



-Trh^{3r-h). 



(51) 



The hi{r) = R — uo{r) is the height of the spherical 
cap C3 and /12(f) = r + £ — uo{r) is the height of the 
spherical cap C2, where uo(r) — {r^ — + R^)/2r is 
the radius corresponding to the plane which contains the 
intersection of the two spheres (Figure [3]) . 

From this the exact solution for the partition function 
and its derivative can be computed to obtain: 

Z{RJ) = ^^2^3^3 _ ^2^2^4 ^ _L^2^6 (53) 



and 



9Z 
OR 



iR,£) 



16 



^2^3 ^2 _ 2Tr^m^ 



(53) 



We shall now consider the more general case of N in- 
teracting solute particles confined in a spherical chamber 
of radius R with a soft- walled potential of the form: 



N 



^(xi,--- ,XAr,i?) = y(xi,--- ,XAr)-f ^$(xfc,i?) 

(56) 



k=l 



where y(xi,-- - jX^r) models the interactions between 
the particles and <i>(xfc, i?) confines the fc*'' solute particle 
to the interior of the chamber, with $ ^ cxd as |xfe| — > R. 

From Eq. and Eq. (HI]), the osmotic pressure can 
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be expressed as: 
NkeT 



P. 



wl 



4^ ^5I^"^^'=^(^l'^2,--.,Xjv) •Xfc)j 

152 (E(-Vx.$(x,,i?)-X, 

\k=l \ 



dR 



where x = (xi,--- ,Xjv) and () denotes the ensemble 
average over the particle configurations weighted by the 
Boltzmann factor exp (— ^'//csT). This was derived by 
making the change of variable x = Rx' in Z and differ- 
entiating in i? = s. 

This expression can be used to characterize the cor- 
rections to the classical theory and the relative contri- 
butions of the microscopic effects of the system. The 
first term is the classical pressure one would expect from 
the van't Hoff's law. The second term arises from the 
particle-particle interactions, and the third term is from 
the particle- wall interactions. It is important to note that 
while we have described the source of each distinct term, 
the terms in fact are coupled by the ensemble average 
which depends on the combination of these effects. 

The second correction term has an intuitive interpre- 
tation as follows. Since the contribution to the pres- 
sure is of the form —S/V ■ x we have that when the par- 
ticle interaction force acting on any particle is toward 
the chamber center, the force acting on the bounding 
wall will be relieved and a negative contribution will be 
made to the pressure. This suggests one mechanism by 
which polymerization reduces osmotic pressure relative 
to free monomers. For a polymer in a typical configura- 
tion, most of the monomers will, for entropic reasons, be 
outside a boundary layer of the wall. As a consequence, 
any individual monomer that makes an excursion toward 
the wall, would on average experience a pulling force to- 
ward the chamber center. From the second correction 
term this will reduce the osmotic pressure relative to free 
monomers. This perspective also provides another way of 
interpreting the wall pressure results for the string model 
from Section IVl 

The third term gives the corrections that arise from 
long-range interactions of the solute particles with the 
boundary wall. In the hard-walled limit this term ap- 
proaches zero. 

We observe that the correction terms to the van't Hoff 
law in Eq. (|57p . other than the partial derivative with 
respect to chamber radius R, take the form of the virial 
from classical mechanics. Virial expansions of the pres- 
sure can be found in several statistical mechanical text- 
books to describe departures of a dilute gas or so- 
lution from an ideal noninteracting particle limit. These 
results are however generally developed within the con- 
text of a large macroscopic chamber, whereas our focus 




4 5 6 7 

L/R = (6k t/k)"^;r 



FIG. 4: Osmotic Pressure of Confined Polymer. Correction 
Factor for a Linear Polymer Confined in a Spherical Cham- 
ber. The classical van't Hoff's Law for A'^ free monomers is 
Pvh = NksT/ ^ttR^ . For stretching an individual bond be- 
tween monomers, we let L — {QknT / KY^'^ denote the length 
at which the bond energy becomes 3k bT. 



is on systems for which interaction length scales become 
comparable to those of the chamber. 



VII. OSMOTIC PRESSURE OF CONFINED 
POLYMERS 

To illustrate how this theory can be applied in prac- 
tice, we now present some numerical results for a model 
of a polymer chain of N monomers confined in a spher- 
ical chamber with a hard-wall confining potential. The 
monomers of the polymer will have coupling given by the 
harmonic bonding energy: 



N 



^(Xl,' 



,Xjv) 



(58) 



J=2 



The osmotic pressure of the system can be computed 
from Eq. (j57p . which in this case reduces to: 



-Pwi — 



NhgT 

1^ 



^(-Vx.F-x,). (59) 



k=l 



We will focus on how the osmotic pressure behaves 
as the bond ing streng th K is varied so that the length- 
scale L — QUbT / K of bond fluctuations varies between 
lengths small and large relative to the chamber size. For 
very small K <C QhsT/R^ the fluctuations in the bond 
length are expected to be very large and the iV monomers 
to behave independently. For very large K ^ QhsT/R^ 
the fluctuations in the bond length are expected to be 
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very small and the N monomers to behave similarly to 
a single particle. Thus in the extreme cases the classical 
osmotic pressures are expected, corresponding to an N 
particle or single particle system. 

To obtain the osmotic wall pressure for intermediate 
values of K, the correction factors were estimated nu- 
merically using the Monte-Carlo method with Metropo- 
lis Sampling [38]; see Figure S) We find that as the 
bond stiffness K ^ 0, the osmotic wall pressure ap- 
proaches the classical van't Hoff osmotic pressure for N 
free monomers. As — > cxd we find that the osmotic 
wall pressure approaches the classical van't Hoff osmotic 
pressure for a single particle. The numerical results show 
how the theoretical framework can be used to capture 
the regime relevant to microscopic chambers, such as a 
polymer confined in a vesicle, in which the bond length is 
comparable to the chamber size. As can be seen from Fig- 
ure|31 the transition between the extremes in the bonding 
strength is gradual and occurs smoothly in K. A simi- 
lar numerical study for the osmotic pressure of polymers 
was conducted within the stochastic immersed boundary 
method in |25,] . 

VIII. CONCLUSION 

We have shown that the osmotic pressure deviates 
markedly from classical macroscopic theories when the 
length-scale of the system becomes sufficiently small so 
that the chamber diameter is comparable in size to the 
length-scale of the interactions of the solute particles. 
In particular, we have explored two ways in which this 
can occur, either through interactions among the solute 
molecules or interactions of the solute molecules with the 
wall. While the classical formulas are not directly appli- 
cable in this regime, we have shown that a theoretical 
framework for the equilibrium osmotic pressures can still 
be developed provided additional microscopic features of 
the system are taken into account. We considered how 



the osmotic pressure is manifested in terms of forces ex- 
erted on the confining wall and in terms of the hydrostatic 
pressure of the fluid in the interior of the chamber. We 
showed through several examples how these notions of 
osmotic pressure can differ for microscopic chambers. In 
particular, we showed how the osmotic wall pressure of 
a two particle string and an N-particle polymer interpo- 
late between the van't Hoff laws associated with the ex- 
tremes of strong and weak bonding, which give effectively 
a single particulate entity or many individual particulate 
entities. 

While the osmotic fluid pressure is related in a simple 
local manner to the osmotic pressure deflned in nonequi- 
librium thermodynamics [23, HI], it behaves a bit less 
intuitively. This was demonstrated for the two particle 
string, where the tight binding limit produced an osmotic 
fluid pressure corresponding to two particulate entities. 
This is in contrast to the osmotic wall pressure which re- 
flected a single confined entity. The distinction appears 
to be related to the local nature of the osmotic fluid pres- 
sure, which can vary throughout the chamber and in par- 
ticular differ between the wall region and center of the 
chamber. Fluid pressure is only expected to be constant 
in mechanical equilibrium when the solute-wall forces in 
the system arc local, but here we have been considering 
interactions between particles and walls on length scales 
comparable to the chamber size. An important conse- 
quence is that the fluid pressure need not be constant 
throughout the chamber interior and the pressure felt by 
the wall need not be equal to the build-up in fluid pres- 
sure in the chamber center. 

The theoretical framework presented here should be 
readily applicable to the study of equilibrium osmotic 
phenomena in many microscopic physical systems aris- 
ing in biology and technological applications. While only 
equilibrium systems where considered here, it may be 
possible to use many of the central ideas of this the- 
ory to ultimately formulate a non-equilibrium (or near- 
equilibrium) theory for microscopic osmotic phenomena. 
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APPENDIX A: OBTAINING THE 
PARTICLE- WALL INTERACTION FORCE AND 
POTENTIAL FROM THE CONFINEMENT 
FORCE AND POTENTIAL (INVERSION 
FORMULAS) 

From Eq. ^ we see that in order to obtain explicit 
expressions for the osmotic pressure it is useful to have an 
expression for G(|x — y|). In modeling systems, we may 
readily have only the confining force F„i(x) while the 
detailed particle- wall interaction forces must somehow be 
inferred. For radial symmetric potentials of the form 
Fwi(x) = i^(|x|)|^ and G(r) = G(|r|)|f|, this requires 
solving the following inverse problem for G given F: 

-F(|x|)^ = / G(|y-x|)^i^dy. (Al) 
|x| Jan |x-y| 

By dotting both sides with this becomes the scalar 
problem: 

-F(|x|) = / G(|y-x|)^:i^.^dy. (A2) 
Jon |x ^ y| 1^1 
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In spherical coordinates, the integral transform can be 

expressed as: 

- F{r) = / G{p){r^ - B? + p')dp (A3) 

where p = |x — y| and r = |x|. 

To ensure a unique solution to the inverse problem we 
shall make the assumption that the particlc-wall inter- 
actions occur only over a distance less than the radius 
R from the chamber wall, that is g{p) = for p > R. 
Under this assumption the transform can be inverted ex- 
actly with the solution formula: 

+p{R-p)F'{R-p)+ [ ' F{s)ds] 



from which G{p) is readily obtained. 

Alternatively, for a given wall potential $(x, i?) the 



inverse problem is to determine a (f){p) so that for all x 
$(x,E) = / </.(|x-y|)dy. (A5) 

In the case that the potential $ is radially symmetric in 
r = |x|, a change of variable allows for the integral to be 
expressed as 

$(r,i?) = — / <i>{p)pdp. (A6) 

JR-r 

A unique solution can be found for this problem under 

the assumption that = for p > i?. By differen- 
tiating both sides in r and substituting p = R — r, the 
following inversion formula is obtained: 

(A7) 



